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I. INTRODUCTION 



The description of doubly heavy baryon properties acquires in the last years the status 
of actual physical problem which can be studied experimentally. The appearance of ex- 
perimental data on B c mesons 0], heavy-light baryons 0] stimulates the investigation of 
heavy quark bound states and can help in discriminating numerous quark models. Recently 
first experimental indications of the existence of doubly charmed baryons were published by 
SELEX |3|. Although these data need further experimental confirmation and clarification it 
manifests that in the near future the mass spectra and decay rates of doubly heavy baryons 
will be measured. This gives additional grounds for the theoretical investigation of the dou- 
bly heavy baryon properties. The success of the heavy quark effective theory (HQET) |3] 
in predicting properties of the heavy-light qQ mesons (B and D) suggests to apply these 



2 



methods to heavy-light baryons, too. The semileptonic decays of heavy hadrons present 
also an important tool for determining the parameters of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. 

The relativistic quark model has been remarkably successful in describing the observed 
hadronic states and their decay rates |E IE 0, 1^ ■ This includes the heavy-light mesons and 
heavy quarkonia. The quark model predicts the existence of doubly heavy baryons containing 
two heavy quarks (cc, be, bb) and one light quark (u, d, s). The energies necessary to produce 
these particles are already reached. The main difficulty remains in their reconstruction since 
these states have in general a large number of decay modes and thus high statistics is required 

Doubly heavy baryons occupy a special position among existing baryons because they can 
be studied in the quark-diquark approximation and the two-particle bound state methods 
can be applied. The two heavy quarks compose in this case a bound diquark system in 
the antitriplet colour state which serves as a localized colour source. The light quark q is 
orbiting around this heavy source at a distance much larger (~ than the source size 

(~ 2 / ttlq) . The estimates of the light quark velocity in these baryons show that its value 
is v/c ~ 0.7 -T- 0.8 and the light quark should be treated fully relativistically. Thus the 
doubly heavy baryons look effectively like a two-body bound system and strongly resemble 
the heavy-light B and D mesons jE E3|- Then the HQET expansion in the inverse heavy 
diquark mass can be performed. We used a similar approach for the calculation of the mass 
spectra of doubly heavy baryons |8j . The ground state baryons with two heavy quarks can be 
composed from a compact doubly heavy diquark of spin or 1 and a light quark. According 
to the Pauli principle the diquarks (bb) or (cc) have the spin 1 whereas diquark (be) can have 
both the spin and 1. 

There exists already a number of papers devoted to studying both the mass spectra of 
doubly heavy baryons and their decay rates 0, IE EH, E3, EE Q EE El, E3, El, EE US HH . 
Whereas the results for the mass spectra are in agreement, the calculations of exclusive 
semileptonic decays lead to essentially different values for the decay rates |IE, EE EE E3 
obtained with the help of the Bethe-Salpeter equation, QCD sum rules and relativistic three 
quark model. 

Here we study semileptonic decay rates of doubly heavy baryons using the relativis- 
tic quark model in the quark - diquark approximation. The covariant expressions for the 
semileptonic decay amplitudes of the baryons with the spin 1/2, 3/2 are obtained in the limit 
m c , nib oo and compared with the predictions of HQET. The calculation of semileptonic 
decays of doubly heavy baryons (bbq) or (beq) to doubly heavy baryons (beq) or (ccq) can be 
divided into two steps (see Fig. Q). The first step is the study of form factors of the weak 
transition between initial and final doubly heavy diquarks. The second one consists in the 
inclusion of the light quark in order to compose a baryon with spin 1/2 or 3/2. 

The paper is organized as follows. In Sec. Owe describe our relativistic quark model and 
present predictions for the masses of ground state heavy diquarks and doubly heavy baryons. 
We apply our model to the investigation of the heavy diquark transition matrix elements in 
Sec. 11111 The transition amplitudes of heavy diquarks are explicitly expressed in a covariant 
form through the overlap integrals of the diquark wave functions. The obtained general 
expressions reproduce in the appropriate limit the predictions of heavy quark symmetry. 
Section IIVI is devoted to the construction of transition matrix elements between doubly 
heavy baryons in the quark-diquark approximation. The corresponding Isgur-Wise function 
is determined. In Sec. IV! semileptonic decay rates of doubly heavy baryons are calculated in 
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Bqq f| pL = Bq/q 

FIG. 1: Weak transition matrix element of the doubly heavy baryon in the quark-diquark approx- 
imation. 



the nonrelativistic limit for heavy quarks. Our conclusions are given in Sec. IVI1 

II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach and quark-diquark picture of doubly heavy baryons the 
interaction of two heavy quarks in a diquark and of the light quark with a heavy diquark in a 
baryon are described by the diquark wave function (^d) of the bound quark-quark state and 
by the baryon wave function (^b) of the bound quark-diquark state, respectively. These 



wave functions satisfy the two-particle quasipotential equation [23, 24] of the Schrodinger 
type 

where the relativistic reduced mass is 

E l E 2 M 4 - (mf - ml) 2 
m ~ E^E 2 ~ UP ' (2) 

and the center of mass energies of particles on the mass shell Ei, E 2 are given by 

M 2 — m\ + m\ M 2 — m\ + m\ 
El = 2M ' E2 = 2M • (3) 

Here M = E\ + E 2 is the bound state mass (diquark or baryon), m 12 are the masses of 
heavy quarks (Q\ and Q 2 ) which form the diquark or of the heavy diquark (d) and light 
quark (q) which form the doubly heavy baryon (B), and p is their relative momentum. In 
the center of mass system the relative momentum squared on mass shell reads 

[M 2 - (mi + m 2 ) 2 ][M 2 - (m, - m 2 ) 2 ] 
1 ] AM 2 ' 1 ' 

The kernel V(p, q; M) in Eq. Q is the quasipotential operator of the quark-quark or 
quark-diquark interaction. It is constructed with the help of the off-mass-shell scattering 
amplitude, projected onto the positive energy states. Here we closely follow the similar 
construction of the quark-antiquark interaction in heavy mesons which were extensively 
studied in our relativistic quark model |26|,|27(. For the quark-quark interaction in a diquark 
we use the relation Vqq = Vqq/2 arising under the assumption about the octet structure 
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TABLE I: Mass spectrum of ground states of doubly heavy baryons (in GeV) |8j. {QQ} denotes 
the diquark in the axial vector state and [QQ] denotes the diquark in the scalar state. 



Baryon 


Quark content 


J p 


Mass 


"cc 


J CC 1/7 


1/2+ 


3.620 


w* 
"CC 


J cc \ct 


3/2+ 


3.727 


a L cc 




1/2+ 


3.778 


i4 CC 




3/2+ 


3.872 




\ uu iH 


1/2+ 


10.202 


"* 


\bb\a 


3/2+ 


10.237 


Qui 

1 L bb 




1/2+ 


1 35Q 


o* 


{66}s 


3/2+ 


10.389 


^cb 


{c% 


1/2+ 


6.933 


"l 

"c6 


[cb]q 


1/2+ 


6.963 


"cfc 


{cb}q 


3/2+ 


6.980 


Qcb 


{cb}s 


1/2+ 


7.088 


Kb 


[cb]s 


1/2+ 


7.116 


Kb 


{cb}s 


3/2+ 


7.130 



of the interaction from the difference of the projection onto QQ and QQ colour states. 
The quasipotential of the quark- ant iquark interaction is the sum of the usual one-gluon 
exchange term and the confining part which is the mixture of long-range vector and scalar 
linear potentials, where the vector confining potential contains the Pauli terms. The explicit 
expressions for these quasipotentials are given in Ref. H. The quark masses have the 
following values m b = 4.88 GeV, m c = 1.55 GeV, m s = 0.50 GeV, m u d = 0.33 GeV. 

We calculated in the framework of the relativistic quark model the mass spectra of heavy 
diquarks and doubly heavy baryon masses in the quark-diquark approximation in Ref. j8j. 
The masses of the ground state axial vector diquarks were found to be M^ v = 3.226 GeV, 
M£ v = 9.778 GeV, M£ v = 6.526 GeV, and the mass of the scalar diquark Mg = 6.519 GeV. 
The calculated masses of the ground state doubly heavy baryons are listed in Table |U 



III. HEAVY DIQUARK TRANSITION FORM FACTORS 

The form factors of the subprocess d{QQ s ) — > d'{Q'Q s )eu where one heavy quark Q s is 
a spectator are determined by the weak decay of the active heavy quark Q — > Q'ev. The 
local effective Hamiltonian is given by 

Heff = ^ QQ , (0'7m(1 - Tb)Q) (Ml - TsK) , (5) 

where Gp is the Fermi constant and Vqq> is the CKM matrix element. In the relativistic 
quark model the transition matrix element between two diquark states is determined by the 
contraction of the wave functions ^ of the initial and final diquarks with the two particle 
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W 



r r« 

FIG. 2: The leading order contribution rW to the diquark vertex function V. 



vertex function T |2 



(d'(Q)\JY\d(P)) 



A, a, p,u = ±- 



1 



(6) 



Here we denote mass, energy and velocity of the initial diquark (QbQ s , index b stands for 
the initial active quark and index s for the spectator) by Mi, Ei = M{V° and v = P/Mi, and 
mass, energy and velocity of the final diquark (Q a Q s , index a means the final active quark) 
by M f ,E f = M f v'° and v' = Q/M f . 

The leading contribution to the vertex function T M comes from the diagram in Fig. |21 
HI HL 111 (we explicitly show spin indices) 



r^(p,q) = r« = tti(q0 7 „(l -75)<(pi)^(q 2 )<(p 2 )(2vr) 3 5(p 2 - q 2 )5° 



(7) 



where the Dirac spinors are 



u\p) 



e(p) + m 



/ 



(^p) 



X 



e{p) 



p2 _|_ m 2^ 



-4 



(8) 



V e (p) + m 

Relativistic four-momenta of the particles in the initial and final states are defined as follows 



pi,2 = tiMv ±Y J n{t \ v )p\ v 



M 



Mi = ei(p) + e 2 (p), 



i=i 

qi,2 = e 1>2 {q)v'±j^n^{v')q i , v' = Mf = ei(g) + e 2 (q) 



i=i 



/ 



(9) 



and are three four vectors defined by 



n®(v) = \v\S t3 + 



:V l V J 



l+v c 

After making necessary transformations, the expression for F should be continued in A4i 
and A4f to the values of initial M$ and final Mj bound state masses. 

The transformation of the bound state wave functions from the rest frame to the moving 
one with four-momenta P and Q is given by 2^, 13 
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*&(<!) = nA^<ex(K) D Vr'( R 0^ ( 10 ) 

where R w is the Wigner rotation, Lp is the Lorentz boost from the diquark rest frame to a 
moving one, and the rotation matrix D 1//2 (i?) is defined by 



'\ 
.0 



JWcw) = 5- i ( Pli2 )5(p)5(p), (ii) 



where 



S ( P ) = J €{p) + m i+ (ap) 



2m y e(p) + m J 

is the usual Lorentz transformation matrix of the four-spinor. 

Using relations JOJ), (fTUj) we can express the matrix element © in the form of the trace 
over spinor indices of both particles. For this aim the following relations are useful 

Sap(A)v%(p) = £ K(Ap)Dl{ 2 (RY p ), 

<T=±l/2 

SJ(P)5S(A) = E ^L 1/2 « P X(Ap). (12) 

<T=±l/2 

Substituting expressions (j2J), (fTUj) in Eq. (JHJ) and using relations (JUJ), (|T2*|) we obtain 

(d'{Q)\Jj\d{P)) = f ^^ (q)^(Q) , lt +m ? , s '^ 

X7m(1 - 7 5 )^p) , =2==^ = «?(0K(0) 
yJ2e b (p)(e b (p) +m b ) 

x y ( ^ + ms) ^(Lq)g(Lp) fe + ms) 



/2e a (q)(e a (q) + m s ) ' y2e a (p)(e s (p) + m a ) 

x^(0)5(p 2 -q 2 )< (p), (13) 

where p b = (e 6 (p),p), p s = (e a (p),-p) and g a = (e a (g),q), q s = (e s (g),-q) (we use the 
notation p = p M 7 M ). The spin wave functions for the scalar and axial vector diquark states 
in the heavy quark rest frame read 0, 0| 

,sr^ 1 + 7° 



n AV (0) = ^7b*, (14) 



where e a is the polarization vector of the axial vector diquark. Using these diquark spin 
wave functions we get the following expression 



(d'{Q)\J™\d{P)) = /^^^,o(q)^o(p)^(P2-q 2 ) 



*r-i, T \ (q s + m s ) T ^s,Av f ^ + m 



xTr S 1 ~\L Q ) ^ ' b> H»> AV (0) 



2e s (q)(e s (q) + m c ) J^a(q)(e a (q) + m a 



7 



x5- 1 (Lq)[ 7m (1- 7 5)] S(Li 



(p b + m b ) 



xH s > AV (0) . { P s + ms)T S T {Li 



^2e b (p)(e b (p) +m b ) 



(15) 



v / 2e s (p)(e s (p) + m s ) 

where ^^(p) is the spin- independent part of the diquark wave function and the superscript 
T denotes transposing. After explicit multiplication of the matrices we get the covariant 
(see Ref. jiij) expression for the transition matrix element 

d 3 p d 3 q r 



(d'(Q)\jW\ d (P)) = 2JM l M_ 



-Tr{V d ,(Q, q) lfl (l - 75 )* d (P,p)}5 d (p 2 - q 2 ), (16) 



where the amplitudes ^ for the scalar (S) and axial vector (AV) diquarks (d) are given by 



e b (p) + m b 



e s (p) + m s 



\ 2e b (p) \ 2e s (p) 
-( 



v + 1 0—1 
+ 



P 2 



2^2 2V2 (e 6 (p) + m b )(e s (p) + m s ) 



v + 1 



V 2V2 e s (p)+m s 



+ 



v-1 



2^/2 e b (p) + m b/ 



V 



7o$s(p), 



(17) 



q AV (P,p,e) 



e b (p) + m b e s (p) + m s 
\ 2e b (p) \ 2e s (p) 
v — 1 2{e ■ p)p 



v + 1 „ v — 1 
■e + 



P 2 



2V2 2V2 (e b (p) + m b ) (e,(p) + m s 
■0 + 1 £p 



2^2 0&0) + m b )(e s (p) + m s ) 2^2 e s (p) + m s 



v — 1 pi 
2\[2 e b (p) + m b 



7o7 5 $ Ay (p). 



(18) 



Here = ^d,o(p) / v2^d is the diquark wave function in the rest frame normalized to 

unity. The four-vector 

( P P)P \ 



p = L P (0,p) 
has the following properties 



(PP) 

M 



P + 



P 



M(E + M) i 



(19) 



( e .p) = -(ep), 

(p-u) = 0. (20) 

The presence of <5 3 (p2 — CI2), with momenta p 2 and q 2 given by Eq. Q, in the decay matrix 
element f)16jl leads to the additional relations 





we s (p) 




w 2 




we s (q) 


% 




w 2 




= we s (p) 




= we s (q) 



iwv'p-Vp) 



Vp 2 



^2 



V WV n -Vim), 



(21) 
(22) 
(23) 
(24) 
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which allow to express either q through p or p through q. The argument of the 5 function 
can then be rewritten as 



e,(p) + £*(?)/ / \ 
p 2 - q 2 = q - P ( v ~ v )' 



(25) 



where w = (v ■ v'). 

Calculating traces in Eq. (JTfij) and using relations one can see that the spectator 

quark contribution factors out in all decay matrix elements. Indeed all transition matrix 
elements have a common factor 



e s (p) +m s 



\ 2e s (p) \ 2e s (q) 



e s (q) +m s 



Iw — 1 



+ 



w + l\e s (p) + m s e s (q)+m, 

VpW 



[e s (q) + m a ][e 8 (p) + m s ] 



w + 1 



I s (p,q). (26) 



If the 5-function is used to express q through p or p through q then I s (p,q) = I s (p) or 
I s (p,q) = T s {q) with 



Up) 



\ 



we s (p) - \Jw 2 - l^/p 2 



e-(p) 



y 1 e s {p)[we s {p) - \/w 2 ^f] 



() \ >Je t (i>) - m 




Iw - 1 
w + 1 



e s (p) + m s 



'w — 1 



w + 1 



e s (p) +m s - J e s (p) - m s . (27) 



The weak current matrix elements have the following covariant decomposition 
(a) Scalar to scalar diquark transition (be — > be) 



(Sfwmstiv)) 



M Si M Sf 



h+(w)(v + t>% + h-(w)(v - u%, 



(28) 



(b) Scalar to axial vector diquark transition (be — > cc 

(AV(v',s'M\S(v)) 



(AV(v>,e')\j£\S(v)) 



h Al (w)(w + l)e'* - h M (w)(v ■ e'*)v^ - h Az (v ■ e'*)v 



VM AV M S Jy ' " 

(c) Axial vector to scalar diquark transition (bb — > 6c) 



(5(TO|J2T|AV(t; |e )> 



ih v (w)e lMX p 1 e a v' p v 1 



(S(v')\j£\AV(v,e)) 



VM AV M S 

= h Al (w)(w + l)e^ - /ia 2 (^)K ■ e)^ - ^a 3 (^ • e') v /o 



(d) Axial vector to axial vector diquark transition (bb — > be, be — > cc 



(29) 
(30) 

(31) 
(32) 



M AVi M AVf 
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+h 4 (w)(v' ■ e)e'; - (i; • e'*)(v' ■ e)[h 5 (w) VlM + foiw)^], (33) 

{AV f {v>,e>))\j£\mv,e)) 



M AVi M AVf 



ie mPl {e^e'^[h 7 {w){v + v') a + h 8 (w)(v - v') a ] 

+v' f> v r [h*(w)(y ■ e'*)e a + h 10 (w)(v' ■ e)e'* a \). (34) 



The transition form factors are expressed through the overlap integrals of the diquark wave 
functions and are given in the Appendix EI These exact expressions for diquark form factors 
were obtained without any assumptions about the spectator and active quark masses. 

If we consider the spectator quark to be light and then take the limit of an infinitely heavy 
active quark mass, m a fc — > oo, we can explicitly obtain heavy quark symmetry relations for 
the decay matrix elements of heavy-light diquarks which are analogous to those of heavy- 
light mesons [lH 3^ 



h + (w) = h v (w) = h Al (w) = h A3 (w) = h Aa (w) = hx(w) = h 3 (w) = h A (w) = h 7 (w) = £(w), 
h-(w) = h A2 (w) = h,A 2 {w) = h 2 (w) = h 5 (w) = h Q (w) = h 8 (w) = hg(w) = h w (w) = 0, 

with the Isgur-Wise function 



««") = V4l / 5) * 3 (p - q + '^'t'-'V - v)) . (35) 

V w + 1 J (zn) 6 \ w + 1 J 

The diquark transition matrix element should be multiplied by a factor 2 if either the initial 
or final diquark is composed of two identical heavy quarks. 

For the heavy diquark system we can now apply the v/c expansion. First we perform the 
integration over q in the form factors (|Al)) - (jA8|) and use relations (|23J). Then, applying the 
nonrelativistic limit, we get the following expressions for the form factors. 

(a) Scalar to scalar diquark transition 

h+(w) = F(w), 

h_(w) = -(w + l)f(w)F(w), (36) 

(b) Scalar to axial vector diquark transition 

h v (w) = [1 + (w + 1) f(w)]F(w), 
h Al (w) = h M (w) = [l + (w- l)f(w)]F(w), 

h A2 (w) = -2f(w)F(w), (37) 

(c) Axial vector to scalar diquark transition 

h v (w) = h A ,(w) = [l + (w + l)f(w)}F(w), 
h Al (w) = [l + (w-l)f(w)]F(w), 

h A2 {w) = 0, (38) 

(d) Axial vector to axial vector diquark transition 

hi(w) = h 7 {w) = F(w), 

h 2 \w) = h 8 (w) = -(w + l)f(w)F(w), 
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1 1.01 1.02 1.03 1.04 1.05 1.06 1.07 
w 

FIG. 3: The function F(w) for the bb — ► be quark transition. 




FIG. 4: The function F(w) for the be — > cc quark transition. 



where 



F(w) 



and 



h 3 (w) = h 4 (w) = (l + {w + l)f(w))F(w), 
h 5 (w) = h 9 (w) = 2f(w)F(w), 
h 6 (w) = h w (w) = 0, 



1/2 



1 + 



W(W+1) \y ^J m 2 + ( w 2 _ X ) m 2 



fW 



(39) 



(27T) 



$F (p+^-(v'-v)) <&,(p) 

V 10 + 1 / 

(40) 
(41) 



The appearance of the terms proportional to the function f(w) is the result of the account of 
the spectator quark recoil. Their contribution is important and distinguishes our approach 
from previous considerations 0, 0|. We plot the function F(w) for bb — > be and be — > cc 
diquark transitions in Figs. lU 
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IV. DOUBLY HEAVY BARYON TRANSITIONS 



The second step in studying weak transitions of doubly heavy baryons is the inclusion of 
the spectator light quark in the consideration. We carry out all further calculations in the 
limit of an infinitely heavy diquark mass, Md — > oo, treating the light quark relativistically. 
The transition matrix element between doubly heavy baryon states in the quark-diquark 
approximation (see Fig. [T] and |2J) is given by [cf. Eqs. (JHJ) and (J7j)] 

(B'(Q)\JW\ B (P)} = 2\ M,.M,. J ^0^*B>, Q (ci)(d'(Q)\ jW\ d (P))* B , P (p)5 3 (P q ~ q,), 

(42) 

where ^ b,p{p) is the doubly heavy baryon wave function; p and q are the relative quark- 
diquark momenta, p q and q g are light quark momenta expressed in the form similar to 
©. 

The baryon ground-state wave function ^b,p{p) is a product of the spin-independent part 
^b(p) satisfying the related quasipotential equation (JTJ and the spin part Ub{v) 

*b,p(p) = *b(p)U b (v). (43) 

The baryon spin wave function is constructed from the Dirac spinor u q (v) of the light spec- 
tator quark and the diquark wave function. The ground state spin 1/2 baryons can contain 
either the scalar or axial vector diquark. The former baryon is denoted by Hqq, and the 
latter one by Hqq/. The ground state spin 3/2 baryon can be formed only from axial vector 
diquark and is denoted by Hqq,. To obtain the corresponding baryon spin-states we use in 
the baryon matrix elements the following replacements 

u a (v) — > Uw> (v), 

[e»(v)u q (v)] spin 3/2 - m* (v), (44) 

where baryon spinor wave functions are normalized by U b IIb = 1 (B = S', S) and the 
Rarita-Schwinger wave functions are normalized by U^iU^n = — 1- 

Then the decay amplitudes of doubly heavy baryons in the infinitely heavy diquark limit 
are given by the following expressions. 

(a) ^'qq s -> S' Q , Qs transition 



(v')\J]?\E' QQs (v)) 



ih + (w)(v + v'), + h_(w)(v-v'),m n (v')U %o (v) V (w), (45) 



( b ) -QQ a E Q'Qs and E QQs -> H' Q , Qs transitions 

- -^[z/iy^Je^^ p w 7 - g tia h A A w + l ) + v^v a h A2 {w) 



2VM 7 M F v^ L 

H^^3H]% Qs K)7 5 (7 a + ^)^ Ws M^(^), (46) 

2JMJ M^ = ^^M^JWr^ w - g m h Al (w + 1) + v^vji^w) 
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+^>A a H]% gs K)(7 a + ^)75f/H ggs (^)^H, (47) 



( c ) E 'qq s ~+ E *Q'Qs and E *QQ S ^ -Q'Qs transitions 



[ih v (w)e fJia p 1 v' p v 1 - g pa h Al (w + 1) + VpV a h A2 (w) 
+v' lx v a h Az {w)\m* qiQs (v')U %Qa (v) V (w), (48) 

[?/iy(w)e MQ/37 t; //3 t; 7 - g m h Al (w + 1) + v^v'Ji A2 {w) 
+v,v' a h A3 (w)}U s , (v')US h (v)r,(w), (49) 



2VM 7 M F 

2^JM I M F 
( d ) -»• S Q'Q S transition 

2y jM lMp = -^{9px[hi{w){v + v) fl + h 2 {w){v-v) ll }-h 3 (w)g fip vx 

-h 4 (w)g llX v' p + v' p v x [h 5 (w)v^ + h 6 (w)v'^] 

+^af3~ / (g p 3 gIMw)(v + v')° + h 8 ( W )(v - v')*] 

+v'^[h 9 (w)gy x + h 10 (w)g a x v']) K)7 5 (7 A + v 



x(Y + vni 5 Uz QQs (vMw), (50) 
(e) S QQs -> S^, Qs and E* QQs -> 2 q ,q s transitions 

-h A (w)g^ x v' p + v'pVxihiw^n + h 6 (w)^] 
+Wt(^aM^ + u') a + M«0(« - v') a ] 
+v'^[h 9 {w)g a p v x + fcioHfl?^])}^^^) 
x( 7 p + ^)75f/H Ws (^H, ' " (51) 



-h A {w)gpxv' p + v'pVxlhiw^f, + h 6 (w)v'p] 
+ie^(g p 3 gl[h 7 (w)(v + v'T + h 8 (w)(v - v'T] 

+v'^[h,{w)g a p vx + h w (w)g%v> p })}u EQlQs (v') 

xi 5 (l X + v' x M h (v) V (w), ' (52) 

( f ) E *QQ S S Q'Q S transition 

{E ^ Q ' { C l MrMp Q ' {V)) = -{^Mw)(v + v')p + h 2 (w)(v-v'),]-h 3 (w)gppvx 

-K{w)gpxv'p + ^aM w )*V + ^e(w)^] 

^(^[^(^(^ + */r + ft 8 w(« - v'T] 

+v > ^[h 9 (w)g«vx + h 10 (w)g a xv'p}))m, Qi Jv') 
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FIG. 5: The Isgur-Wise function rj(w) of the light quark - heavy diquark bound system. 



xU£. (v)r](w). 



(53) 



Here rj(w) is the heavy diquark - light quark Isgur-Wise function which is determined by 
the dynamics of the light spectator quark q and is calculated similarly to Eq. (|33j) 



r]{w) 



w + 1 



(v 



(54) 



We plot the Isgur-Wise function r](w) in Fig. El In the nonrelativistic limit for heavy quarks 
the diquark form factors hi(w) obey relations (|36j) - (j39j) . In this limit the baryon transition 
matrix elements contain the common factor F(w)r)(w) (cf. 0]). 



V. SEMILEPTONIC DECAY RATES OF DOUBLY HEAVY BARYONS 



The exclusive differential rate of the doubly heavy baryon semileptonic decay B — > B'ev 
can be written in the form 



dT Gl\V bc \ 2 M 3 F Vw 2 - 1 



dw 



48tt 3 



(55) 



where w = (v ■ v') = (Mj + M| - £; 2 )/(2M 7 M F ), k = P - Q and the function Q(w) is the 
contraction of the hadronic transition matrix elements and the leptonic tensor. For the 
massless leptons the differential decay rates of the transitions z!qq s —> ^q'q s and 5qq s — ► 
Sg,g s are as follows 



dw [ ^ QQ ° 



J = ^ifif'V - 1 



1/2. 



X 



72tt 3 

, 9 / \ W ~ 1 7 • 



tw + 1) 3 M]U 2(M 7 2 + Mj. - 2M I M F w) 



w + 1 



h 2 v H 



+ 



(M F - M T w) 



x/i Al (iw) + («;-!) (M F h A2 (w) + Mjh A3 (w)) 



>T] 2 (w), 



(56) 



Mfd^^ QQs 



t ,Mj^M F ,h A2 



hA 2 ) 



(57) 
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TABLE II: Semileptonic decay rates of doubly heavy baryons 5^ and (4 n x 10 14 Gev). 



Decay 


our 


Ref.fl6] 


Ref.[15j Ref.[17] 


Ref.fl3] 


-bb - 


"* "6c 


1.64 


4.28 






"bb - 


■* "6c 


3.26 


28.5 


8.99 




"bb - 


~* "6c 


1.05 


27.2 


2.70 




"bb 


"6c 


1.63 


8.57 






"* 
"bb 


1 — oc 


0.55 


52.0 






"bb 


"* "6c 


3.83 


12.9 






"be ~ 


"CC 


1.76 


7.76 






"be ~ 


"CC 


3.40 


28.8 






-be - 


"CC 


4.59 


8.93 


4.0 8.87 


0.8 


^bc - 


"CC 


1.43 


14.1 


1.2 2.66 




"be 


"CC 


0.75 


27.5 






"* 

"be 


"CC 


5.37 


17.2 






TABLE III: Semileptonic 


decay rates of doubly heavy baryons il^ and 0,f, c (in 


xl0~ 14 Gev). 


Decay 




r 


Decay 


r 


^66 - 


^'bc 




1.66 


O' —> O 

a L^ c r i L cc 


1.90 


^66 " 


-> ribc 




3.40 


O' — v o* 

"be ^ "cc 


3.66 


^66 " 


-^bc 




1.10 


^bc > ^cc 


4.95 


O* - 


^'bc 




1.70 


^bc — >■ ^cc 


1.48 


O* - 
Si 66 






0.57 


O* — > o 


0.80 


O* - 
Si 66 


^ "be 




3.99 


"be ^ "cc 


5.76 



The differential decay rates for other transitions are given in the Appendix iBl 

The semileptonic decay rates of doubly heavy baryons are calculated in the nonrelativistic 
limit for heavy quarks and presented in Tables HI and IIHI 

VI. CONCLUSIONS 

In this paper we calculated the semileptonic decay rates of doubly heavy baryons in the 
quark-diquark approximation. The weak transition matrix elements between heavy diquark 
states were calculated with the self-consistent account of the spectator quark recoil. It 
was shown that recoil effects lead to the additional contributions to the transition matrix 
elements. Such terms were missed in the previous quark model calculations. If we neglect 
these recoil contributions, the previously obtained expressions [3, HH for heavy diquark 
transition matrix elements are reproduced. It was found that these recoil terms which are 
proportional to the ratio of the heavy spectator to the final active quark mass [see Eqs. ()36|) - 
(jS| ] give important contributions to transition matrix elements of doubly heavy diquarks 
even in the nonrelativistic limit. In this limit these weak transition matrix elements are 
proportional to the function F(w) (|4Uj) which is expressed through the overlap integral of 
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the heavy diquark wave functions. The function F(w) falls off rather rapidly, especially for 
bb — > be diquark transition where the spectator quark is the b quark (see Figs EI). Such a 
decrease is the consequence of the large mass of the spectator quark and high recoil momenta 
(<?max ~ m b — m c ~ 3.33 GeV) transfered. 

We calculated the doubly heavy baryon transition matrix elements in the heavy diquark 
limit. The expressions for transition amplitudes and decay rates were obtained for the most 
general parameterization of the diquark transition matrix elements. The Isgur-Wise function 
rj(w) (J51j) for the light quark - heavy diquark bound system was determined. This function 
is very similar to the Isgur-Wise function of the heavy-light meson in our model [28| as it 
is required by the heavy quark symmetry. In the heavy quark limit the baryon transition 
matrix elements contain the common factor which is the product of the diquark form factor 
F(w) and the Isgur-Wise function 77(10). 

Our results for the semileptonic decay rates of doubly heavy baryons E bb and E bc are 
compared with previous predictions in Table It is seen from this Table that results of 
different approaches differ substantially. Most of previous papers give their 

predictions only for selected decay modes. Their values agree with our in the order of 
magnitude. Our predictions are smaller than the QCD sum rule results ^3] by a factor of 
~ 2. This can be a result of our treatment of the heavy spectator quark recoil in the heavy 
diquark. On the other hand the authors of Ref. ^6|, where the Bethe-Salpeter equation is 
used, give more decay channels. Their results are substantially higher than ours, for some 
decays the difference reaches almost two orders of magnitude which seems quite strange. 
E.g., for the sum of the semileptonic decays E bb — > Ref. [l6j predicts ~ 6 x 10~ 13 which 
almost saturates the estimate of the total decay rate r|°^ tal ~ (8.3 ± 0.3) x 1CT 13 [2] and 
thus is unlikely. 

All considerations in the present paper were done in the heavy quark limit. The calcu- 
lations of the semileptonic decay rates of heavy mesons indicate that the 1/rriQ corrections 
to the heavy quark limit (especially l/m c ) are important Thus they can give 



sizeable contributions also to the semileptonic decay rates of doubly heavy baryons. Their 
account will result in the more complicated structure of the baryon matrix elements. We 
plan to study these corrections in future. 
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APPENDIX A: FORM FACTORS OF THE DIQUARK TRANSITIONS 



h+(w) = h(w) = h 7 (w) = */_=- / ^^Mq 



w + lJ (2tt)3 2e a (q) \ 2e b (p) 



e a (q) + m a 



e b (p) + m b 
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x 1 



[e a (q) +m a ][e h (p) +m b ] 



MP) 



xl s (p,q) 8* p- q + 



w + 1 



(v' - v) , (Al) 



= h 2 (w) = h$(w) 



d 3 p d 3 q - , 



w + l J (2tt) 3 2e a (g) \| 2e 6 (p) 



e a (g) + m a 



e 6 (p) + m b 



x> 



lw + 1 

w - 1 \ e a(<?) + m a £b(p) + rn b/ 



Mp) 



XI sip, q) S 3 f p - q + - v) ) , (A2) 



h v (w) = h 3 (w) = h 4 (w) 



2 /■ d 3 p d 3 g - . 



x 



tu + 1 
1 + 1 



(2tt) 3 

W + l f \^q I 



e a (g) +m a 



e 6 (p) + m fe 



^ 2e a (g) \ 2e b (p) 



+ 



w - 1 \e a (q) + m a e b (p) + m b J 



[e a (q) + m a \[e b {p) + m b ] 
xl s (p,q) S 3 (p- q + 



W + l 



V — V 



(A3) 



/lAi(w) = 



d 3 p d 3 q 
w + lJ (2tt) 3 



x 



1 + 



'u> — 1 



Mq) 



e a (?) + ™ a 



e 6 (p) + m 6 



^ 2e a (g) ^ 2e 5 (p) 



+ 



w + 1 \ e a(g) + m a e b (p) +m b J [e a (q) + m a \ [e b (p) + m b \ 



x<Z> d (p)I s (p,q) 5 3 p-q + 



3/„ „ , e s(p) + £*{<!) ( ,j 



w + l 



v'-v) , 



(A4) 



h A2 {w) = h 5 (w) = h 9 (w) = 



x ■ 



d 3 p d 3 q 
w + lJ (2tt) 3 

2v^ 



e a (q) + m a 



\ 2e a {q) \ 2e b {p) 



e b (p) + m b 



iw 



i Vp 2 



\lw 2 - l[e a (q) + m, 
Xl s (p,q) Wp-q + 



:t/ , € s(p) + e s (q),_j 



w + 1 



w-1 e b (p) + m b 
v' " V)) , 



*d(p) 



(A5) 



^(w) = h 6 (w) = h w (w) 



d 3 p d 3 q - 



W + iy (2tt) 3 dKHJ \ 2e a (q) \ 2e b (p) 



e a (q) + m a 



e b (p) +m b 
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x 



'w + l 



Vw 2 - l[e b (p) + m b ] 



w- le a (q) +m a 



xl s (p, q) 5 3 p - q + ^ -_^ ( v ' - v) ] , (A6) 
\ w + l 



_ /■ d 3 pd 3 q - e a (q) + m a e b (p) +1 

lJ (2ir) 3 d ' lq ^ 2e a (q) \ 2e b (p) 

/w + l (w — I y/q 2 y/p 2 

V w - 1 \w + 1 e a (g) + m a e b (p)+m b/ 

w + 1 V^Vq 1 



w + 



e 6 (p) + m b 



x 



1 + 



w - 1 [e a (q) + m a ] [e b (p) + m 6 ] 



*d(p) 



x/ a (p,g) 5 3 p 



e,(p) + e( g ) 

w + l v 7 , 



(A7) 



^ 3 (w) 



w + l 



d 3 pd 3 q ^ 

W d ' (q) ' 



-o(5) 



e 6 (p) +m b 



x 



'w + l 



1 + ., . 

V w - 1 \e a {q) +m a 

w + l Vj^Vq 2 



+ 



w - 1 [e a (g) + m a }[e b (p) + m b ]_ 
p-q + 



\J 2e a (q) \ 2e b (p) 
w — 1 \ 
w+ le 6 (p) +m b J 

*d(p) 



w + l 



V — V 



(A8) 



APPENDIX B: DIFFERENTIAL DECAY RATES OF DOUBLY HEAVY 

BARYONS 



dw [ ^ QQ ° 



9MY99f (w 2 _ 1)3/2(w + 1)M 3 (M/ + Mf)2 

247T ;i 

Mi- M F 



X 



/i+(w) 



M 7 + M F 



77 (w), 



(Bl) 



dw 



^|VQQf (w2 _ 1)1/2(w )3M 3 

36tt 3 V ; V ; f 



x <^ (M 7 2 + M| - 2M I M F w) 



Sh 2 Al (w) 



w 



1 



w + l 



+ (w- 1) 



M 7 2 (w + l)h 2 Al (w) + (w - 1) {M F h M {w) + M / /i As (w)) 5 



+2{M F - M lW )h Al {w) (M F h A2 (w) + Mj/UsH) 



>77 2 (w), (B2) 
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1 M 3 dT ~ 
Q'Qs) = 2 J^dw^'^ ~* E h'Qs> M i ^ M F> h A 2 -»• h A2 ,h A3 -> /i A3 ), (B3) 



x (j2 - w) \±h 2 {w) -(w- l)(hg(w) + h w (w)) 
vo — X 

+^y[/isH + h 4 (w)f + (w- l)(2h 7 (w) + hg(w) + h 10 (w))'- 

W — 1 . . , . , ,l2 



+4 

ID 



(Mj - M F )h 7 (w) 
- 1 



iu + 1 



(Mj + M F )h 8 (w) 



(w + 2)[(M 7 + M^h^w) - (M 7 - M F )/i 2 H] 



4 w + 1 

+ (M/ - M F w)h 3 (w) + (M F - Miw)h 4 (w) 
+ {w 2 - l)[M F h 5 (w) + Mjheiw)}) W(w), 



(B4) 



) = 



1087T 3 



l) 1/2 (u> + 1) 3 M£<| (Mj? + M| — 2MiM F w) 



\3 !i /r3 



x /i 2 (w;) + [/i 7 (w) + (w - l)(/i 9 H + /iio(w)] : 



u; - 1 



(/i 2 (w) + [/i 8 (iu) + (u> + l)(/i 9 H - M™)] 2 ) 



+ 3 , 

\w + 1, 

+2^y([/i 3 (w) - h 4 (w)} 2 + h 3 (w)h 4 (w) 
-{w 2 - l){[h 9 {w) - h 10 {w)] 2 + h 9 (w)h 10 (w)) 

yj J. -| 2 

(Mi - M F )h 7 (w) -(Mi + M F )h 8 (w) 

w + 1 J 

^^^( (w - l)[(Mj + M f )/hH - (Mj - M F )/i 2 H] 



+ 



u; + 1 

+(Mi - M F w)h 3 (w) + (M F - Miw)h 4 (w) 
+ (w 2 - l)[M F h 5 (w) + MMw^YWiw), 



(B5) 



i rfr ,„ 

2dw 



y-QQs 



(B6) 



w 2 - iy /2 (w + If M£ j(M/ + - 2MiM F w) 

2\ 



G 2 f \Vqq'\' 

2167T 3 

x (^5(/i 2 7 H + [/i 7 H + (w- l)(h 9 (w) + h 10 (w)} 2 ) 

+3 (^j) 2 (^sH + [fcsH + (w + i)(V«0 - M«0] 2 ) 

+2^— ^2[/i 3 (w) - h 4 {w)f + 5h 3 {w)h 4 {w) 
w + 1 v 
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-(w 2 - l)(2[h 9 (w) - h w (w)] 2 + Bh 9 (w)h 10 (w)))^ 

+5[(M/ - M F )h 7 (w) - ^^y(M/ + M F )h s (w) 

+^T[(( 2w + l )[( M i + M F )hx{w) - (M T - M F )h 2 (w)} 

+ (M 7 - M F w)h 3 (w) + (M F - Mjw)h4,(w) 

+ (w 2 - l)[M F h 5 (w) + Mih 6 (w)]\ + ((Mi - M F w)h 3 (w) 

\ 2 

+ (M F - Miw)h 4 (w) + (w 2 - l)[M F h 5 (w) + Mih 6 (w)]\ 
-2(w + 2)(w- 1) ((Mi + M F )h 1 (w) 

-(M I -M F )h 2 (w)^ V 2 (w). (B7) 
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